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NOTATION 


_  ■  *  .i 


Symbol 

Dimensionless  Form 

Definition 

B 

o.  B 

Buoyancy  force,  positive  upward 

ca 

Center  of  buoyancy  of  submarine 

cc 

Center  of  mass  of  submarine’ 

i 

X 

!, 

*  '  « -r—i 

*  K 

Moment  of  inertia  of  submarine  about  x  axis 

1 

y 

ly 

V  *  Rr 

Moment  of  (inertia  of  submarine  about  y  aids 

i 

i 

1 

r  ihi* 

Moment  of  inertia  of  submarine  about  r.  axis 

i 

*y 

I  1 » 

*y 

Product  of  inertia  about  xy  axis 

i 

y* 

i  1  « lx* 
y*  iuT 

Product  of  inertia  about  yx  axes 

I 

IK 

.  '** 

*« ' 

Product  of  inertia  about  cx  axes 

K 

K'.*, 

to<  U1 

Hydrodynamic  moment  component  about  x 
axis  {rolling  moment ) 

K. 

K. 

K- 

Rollins  moment  when  body  angle  {*,  5)  and 
control  surface  angles  are  aero 

% 

K  •  ■  Kt»  ■ 

•n  i^Tu1 

Coefficient  used  it?  representing  K*  as  a 
function  of  (tf- 15 

K 

P 

K  '  ■  - 
P  $pt4U 

Firet  order  coefficient  used  la  representing 
Keen  function  of  p 

k.*.JSL 

M* 

Coefficient  used  la  representing  K  a#  a  function 
of  $ 

Kplpl 

kpipi‘  ■ 

Second  order  coefficient  used  in  representing 

K  es  s  function  of  p 

k« 

Coefficient  need  in  representing  K  as  a  function 
of  the  product  pq 

ti 


yn 


K 

'I*' 


K' 


T 


K 


v<t 


wp 


“6r 


U 


M 


P 9 


M 

9 


U. 


u. 


r.  '  «  — Sf. 

*r  W 

Kr  '  JOT 

Kr  X$T 

„  i  *v 

Kv  ‘JOT 

K* 


K. 


K  '  «  , 

V<1  *(Vt 


kot4 

KVn 


pi* 
UP  *pl‘ 


K~, 

K  '  « 


**  *ol4 

KSr 

6r  $p?  I 7 

V  »  J 


4ot' 


U 

JOT* 

u* 

•JOT 

U  ' 

pp 

if 

»  ^ 

M 

d 

it>t4U 

u 

<5*1 

’trft 

Mr' 

d 

•1^ 

Coefficient  used  in  representing  Kate 
function  of  the  product  qr 

Tint  order  coefficient  used  in  representing 
K  t*  a  function  of  r 

Coefficient  used  in  representing  K  at  a 
function  of  {• 

Tirat  order  coefficient  used  in  representing 
Kata  function  of  v 

Coefficient  used  in  representing  K  as  • 
function  of  i 

Second  order  coefficient  used  in  representing 
K  as  a  function  of  v 

Coefficient  used  in  representing  Kali  function 
of  the  product  vq 

Coeffic.ent  used  in  representing  K  as  *  function 
of  the  product  vw 

Coefficient  used  is  representing  K  as  a  function 
of  the  product  wp 

Coefficient  used  is  representing  K  as  a  function 
of  the  product  *r 


First  order  coefficient  need  is  representing 
K  as  a  function  of  tf 

Overall  length  of  subcuartaa 

his  as  of  subasarius,  including  water  is  fr as  - 

flooding  spaces 

Hydrodysumic  asomsst  compumt  shout  y  sals 

(pttchiug  momma*) 

Pitching  tnceucus  when  body  angles  («,  9)  sod 
control  surface  angles  are  aero 

Sec oW;  order  coefficient  jssd  is  representing 
Mesa  function  of  p.  First  order  coefficient  is 
ssro. 

First  order  coefficient  used  is  reprase  ruing 
hi  as  a  function  of  q 

First  order  coefilcisur  used  iu  ropresesstug 
as  a  function  of  fq-  I ) 

Coefficient  «iaf  is  vepreseetisg  has* 
function  of  q 


“dhl 

“<»hl 

Second  order  coefficient  used  in  representing 

U  !•  i  function  of  q 

Mhl 

■ ,  J^JaLL* 

*•  *pt*u 

Coefficient  uecd  In  representing  ss  s 

function  q 

% 

M  ' 

n> 

fpZ^ 

Coefficient  used  in  representing  M  ss  a 
function  of  the  product  rp 

U 

rr 

M 

''P 

M  ' 
rr 

M  1 
VP 

Mrr 

*  -s^r 

Mvp 

Hoi* 

Second  order  coefficient  used  in  representing 

M  ss  a  function  of  r.  First  order  confident 
is  tero 

Coefficient  used  in  representing  Mats 
function  of  the  product  vp 

M 

vr 

M  ' 

vr 

Mvr 

$0t4 

Coefficient  used  in  representing  Mm 
function  o‘  the  product  vr 

M 

vv 

M  ' 
vv 

Mvv 

’*^7T 

Second  order  coefficient  used  in  representing 

M  ss  a  function  of  v 

M 

w 

M  ’ 

» 

First  order  coefficient  used  in  representing 

M  ss  a  function  of  w 

F&t-i; 

M 

wtj 

M 

WV) 

tot  v 

First  order  coefficient  used  In  representing 
ss  a  'unction  of  (*t-l) 

%A. 

Coefficient  used  in  representing  Mass  function 
of  « 

w 

'  n^- 

M|* 

,  ^1  w  1 

I 

First  order  coefficient  used  in  representing  M 
ss  s  function  of  w;  equal  to  sero  for  symmetries! 
function 

M|wh 

MU  i 

- 1  Ml*  <<» 

<i  F^ 

Coefficient  used  in  representing  Mq  ss  s  function 
of  « 

u  1  1 

.*  .  _  M»  i*  1 

Second  order  coefficient  used  in  representing 
Mass  function  of  w 

First  order  coefficient  used  in  representing 
**w]  w  |  B*  a  functiott  (t»*D 

M 

wnm 

U 

'Tu^ 

Second  order  coefficient  used  lo  representing 
Mass  function  of  w;  equal  to  sere  for  oyro- 
avetrlcsl  function 

U6b 

M6t' 

M6b 

jpc^u1 

First  order  cci '.f’.cient  used  in  representing 
Mass  function  of  6^ 

Mdo 

..  “5. 

'mTu7 

First  order  coefficient  used  in  representing 
Um  function  of  ft  u 

u6.n 

^8»n 

First  order  coefficient  used  lo  representing 

M|c  sea  .'unction  of  (*>-!} 

Iv 


N. 


N 


N, 

P 


M 

pq 


q*- 


N 


N'  ’ipPT" 
N* 

N*‘ “T^rni* 

N„ 


N  ' 


p  *  $p<«  U 

N.'  *  T-^T- 
P 


tq 


Nr 


pq  4c>4& 

Nqr 

N  '  e  ,-0- 

q*-  $p£s 


N 

N  '*  - I— 

T  £pi*U 


rn 


n  ■ .  Sn l 

rrj  tp-(«U 


N. 

r 


N- 

M.1  *  - — *— 
r 


n 


N 


N 


rlrl 


N, 


N  •  e 
v 


"rlrl 
,  ^  Nl  r  I  6i 

*•> 

N. 


N  ■  *  — 2L 
V’T  fct'U 


N, 


N.'  - 


bi* 


N 

N  1  .  -1v3.. 

'•q  M4 

X,  , 

'  v '  r 

w  ,  _  Nl »!r 

Irlr  * 

N  ,  . 

..  .  Krlr! 

Nvlv| 

v  1  v  1 

Nr]v|, 

N  i  r 
v  (  v  1  q 

Hi  |  1  ii  * 

vM*  foe 

Hydrodynamic  component  about  t 

axis  (yawing  momen'l 

Yawing  moment  when  body  angle*  («,  0)  and 
c  tttrol  iurlace  angle*  are  aero 

Fir*t  order  coefficient  used  In  representing  N 
a*  a  (unction  o(  p 

Coefficient  u*ed  in  representing  N  a*  a  (auction 
of  p 

Coefficient  u*ed  in  representing  N  a*  a  function 
of  tiie  product  pq 

Coefficient  u*ed  in  representing  N  a*  a  function 
of  the  ps  oduct  qr 

Firat  order  coefficient  used  in  representing  N 
a*  a  function  of  r 

First  order  coefficient  used  in  representing 
Nr  as  a  function  of  (tj-1) 

Coefficient  used  in  representing  N  as  a  function 
of  r 

Second  order  coefficient  used  in  representing 
N  as  a  function  of  x 

Coefficient  used  in  representing  N.  at  a 
function  of  r 

First  order  coefficient  usod  in  representing  N 
as  a  function  a  v 

First  order  coefficient  ueed  in  representing  H*. 
as  a  function  of  (tv  1 ) 

Coefficient  used  in  repreeeating  N  as  a 
function  of  v 

Coefficient  used  in  representing  Nssa  function 
of  the  product  vq 

Coefficient  ueed  in  representing  Nr  as  a 
function  of  v 

Second  order  coefficient  need  ia  rsq  resenting 
N  as  e  function  of  v 

Firet  order  coiffieieet  used  ia  fepre*««xiag 
Ny|  vj  ««  a  function  of  (q*  11 

e 


Coefficient  used  in  representing  N  as  a  function 
of  the  product  vw 

Coefficient  used  in  representing  N  as  a  function 
of  the  product  wp 

Coefficient  used  in  representing  N  ss  a  function 
of  the  product  wr 

Firnt  order  coefficient  used  in  representing  N 
as  s  function  of  6r 

First  order  coefficient  used  in  representing 
N6r  ••  »  function  of  (tj*  I) 

Angular  velocity  component  About  y  Axis 
relative  to  fluid  (roll) 


Angular  Acceleration  component  About  x  axis 
relative  to  fluid 


Angulxr  velocity  component  about  y  axis  relative 
to  fluid  (pitch) 


Angular  acceler  -tlon  component  about  y  axis 
relative  to  fluid 

Angular  velocity  component  About  *  Axis 
relative  to  fluid  (yaw) 


Angular  Acceleration  component  about  z.  axis 
relative  to  fluid 


Linear  velocity  of  origin  of  body  axes  relative 
to  fluid 

Component  of  U  in  direction  of  the  x  axis 


Time  rate  of  change  of  u  in  direction  of  the 
x  axis 

Command  speed:  steady  value  of  ahead  epeed 
component  u  for  a  given  propeller  rpm  v  hen 
body  anglee  (a,  d)  a od  control  surface  angles 
are  sero,  Sign  changes  with  propeller  reversal 

Componen*  of  U  in  direction  of  the  y  axle 


Time  rat*  of  change  of  v  in  direction  of  the 
y  axis 

vi 


r 


/ 


W 

w  s  — 

u 

Component  of  U  in  direction  of  the  a  exit 

w 

:  ,  *4 

W  =  - - 

u1 

Time  r»tr  of  change  of  w  in  direction  of  the 
i  axis 

W 

*  =^PiF 

V'eight,  including  water  in  free  flooding  apace* 

X 

X 

X  ®  - 

4 

Longitudinal  body  axle;  alao  the  coordinate  of  a 
point  relative  to  the  origin  of  bed)  uxea 

*B 

v-? 

The  x  coordinate  of  CB 

The  x  coordinate  of  CC 

*0 

,  *o 

*°  s  T 

A  coordinate  of  the  displacement  of  CO  relative 
to  the  origin  of  a  aet  of  fixed  axe* 

X 

.  X 

"ToPt? 

Hydrodynamic  fotce  component  along  x  axi* 
(longitudinal,  or  axial,  force) 

X 

qq 

X  '=3sL 
qq  $p/.4 

Second  order  coefficient  uaed  in  representing 

X  at  a  function  of  q.  First  order  coefficient 
ia  aero 

X 

ro 

x  '  =  —IS. 
rP  if** 

Coefficient  u*ed  in  representing  X  aa  a  function 
of  the  product  rp 

X 

i  r 

X. 

u 

y  .  *rr 

Xrr  St?r 

xu  -isr 

Second  order  coefficient  used  in  representing 

X  a*  •  function  of  r.  Firat  order  coefficient  i* 
aero 

Coefficient  used  in  representing  X  aa  a  function 
of  u 

X 

uu 

X  '  a  ^ 
uu  tc* 

Second  order  coefficient  used  in  repreeeating 

X  aa  a  function  of  u  in  the  non* propelled  cane. 
First  order  coefficient  is  aero 

X 

VT 

xvr 

x  ■  »  T-^ 

VT  *£*’ 

Coefficient  used  in  representing  X  a*  •  function 
of  the  product  vr 

X 

w 

xvv 

Second  order  coeffient  used  in  representing  X 
as  a  function  of  »,  Piral  order  coefficient  ia  aero 

^e*5 

First  order  coefficient  used  ia  representing 
ene  function  o»'  fiy-1) 

X 

•q 

X 

«q  tp4* 

Coefficient  need  la  representing  X  as  a  function 
of  the  product  eq 

vii 
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wW 

X~w 
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u 
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Xib6b 
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V  .  _  x&r$ T 

°r4r 
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V* 

? 

y  •  ^ 
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Y. 
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Y 

P 

TP 

rp 

Y, 

P 

Y.- « JL 

P  Itol* 

Vp‘ 

Y  .  •  » 

p!p!  tpf,* 

Second  order  coefficient  ueed  in  representing 
X  as  »  (unction  of  w.  Fire*  order  coefficient  is 
zero 

First  V'^r  coefficient  used  in  representing  Xy^ 
is  s  function  of  (r? - 1 ) 

Second  order  coefficient  used  in  representing  X 
ss  s  function  of  6^.  First  order  coefficient 
is  zero 

Second  order  coefficient  used  in  representing 

X  as  s  function  of  6r.  First  order  coefficient  is 
zero 

First  order  coefficient  used  in  representing 
Xjr£r  4*  4  function  of  (tj-l) 

Second  order  coefficient  used  in  representing  X 
*•  s  function  of  First  order  coefficient  is 
zero 

First  order  coefficient  used  in  representing 
X,  .  at  function  of 

Latere!  body  axis;  also  the  coordinate  of  a 
point  relative  to  the  origin  of  body  axes 

The  y  coordinate  of  CB 


The  y  coordinate  of  CG 

A  coordinate  of  the  displacement  of  CG  relative 
to  Che  origin  of  a  set  of  fixed  axes 

Hydrodynamic  force  component  along  y  axis 
(lateral  fores) 
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First  order  coefficient  ueed  in  representing  Y 
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ABSTRACT 


''Standard  equations  of  motion  are  presented  for  use  in  submarine 
simulation  studies  being  conducted  for  the  U.S.  Navy.  The  equations  are 
general  enough  to  simulate  the  trajectories  and  responses  of  submarines 
in  nix  degrees  of  freedom  resulting  from  various  types  of  normal 
maneuvers  as  well  as  for  extreme  maneuvers  such  as  those  associated 
with  emergency  recoveries  from  mternplane  jam  and  flooding  casualties. 
Information  is  also  presented  pertaining  to  the  hydrodynamic  coefficients 
and  other  input  data  needed  to  perform  simulation  studies  of  specific 
submarine  designs  with  the  Standard  Equations  of  Motion. 
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The  development  of  the  NSRDC  Standard  Equations  of  Motion  for 
Submarines  and  associated  experimental  techniques  to  provide  the 
necessary  hydrodynamic  coefficients  was  sponsored  primarily  by  the 
General  Hydromechanics  Research  Program  (Sub  Project  SR  009  01  01 
Task  0102)  Additional  support  was  provided  by  the  Submarine  Safety 
Program  (Sub  Project  S4611010  Tasks  11077  and  11081). 


INTRODUCTION 

The  Naval  Ship  Research  and  Development  Center  was  requested  to 
ovide  the  Naval  Ship  Engineering  Center  with  a  report  on  the  general 
equations  of  motion  currently  being  used  for  submarine  simulation  stud¬ 
ies.1  The  primary  purpose  of  the  report  is  to  establish  standards,  pro¬ 
vide  guidance,  and  establish  a  firm  basis  for  contract  negotiations.  The 
report  would  thus  serve  to  facilitate  the  Navy's  dealings  with  contractors 
engaged  in  the  manufacture  of  equipment  such  as  training  simulators  and 
automatic  control  systems  for  submarines.  Particular  interest  was 
expressed  in  having  equations  which  would  fully  describe  the  hydrody¬ 
namic  aspects  of  a  submarine  experiencing  a  casualty  while  In  ahead 
motion,  with  provisions  for  inserting  detailed  representations  of  the 
casualties  and  recovery  systems. 

The  Stability  and  Control  Division  of  the  Hydromechanics  Laboratory 
has  been  conducting  "submarine  motion  simulation  studies  tor  *  number 
of  years  in  connection  with  its  various  assigned  responsibilities.  Such 
studies  have  been  utilised  with  considerable  success  to  solve  a  wide 
variety  of  problems  pertaining  to  the  design  and  operation  of  submarines 
from  the  standpoint  of  stability  end  control.  Most  recently,  these 
studies  have  included;  emergency  recovery  from  sternplane  jams,  loao 
supportability,  and  vertical  ascents  after  emergency  blowing  (both  with 
and  without  forward  propulsion).  These  simulation  studies  have  been 
made  possible  as  the  result  of  a  progressive  in-house  development  of 


i  References  are  listed  on  page  27  . 


six-degree-of -freedom  motion  equations  and  a  coordinated  development  of 
advanced  experimental  and  analytical  techniques  to  provide  the  necessary 
coefficients. 

The  equations  of  motion  as  they  now  exist  at  NSRDC  are  general 
enougn  to  simulate  trajectories  of  submarines  for  all  of  the  various  types 
of  normal  maneuvers  and  emergency  situations  that  have  been  of  concern 
up  to  the  present  time.  Furthermore,  it  is  believed  that  they  will 
adequately  cover  most  of  the  situations  that  may  be  reasonably  anticipated 
for  the  future.  These  equations,  which  form  the  subject  of  this  report, 
have  been  used  continuously  for  the  past  three  years  with  little  change, 
and  are  considered  to  be  standard  at  NSRDC.  In  the  interest  of  consistency, 
it  is  proposed  herein  that  these  equations  be  adopted  and  serve  aa  the 
U.S.  Navy  Standard  until  notice  is  given  by  NSRDC  that  they  should  be 
changed.  It  should  be  emphasised,  however,  that  theae  equations  have 
been  validated  only  with  the  types  of  coefficients  being  generated  by  the 
Hydromechanics  Laboratory.  Therefore,  prospective  users  are 
advised  to  first  contact  NSRDC  to  determine  whether  the  required 
hydrodynamic  coefficients  are  available. 

It  is  the  intention  eventually  to  issue  a  complete  and  comprehensive 
report  on  the  NSRDC  Standard  Equations  of  Motion  for  Submarines  which 
would  be  tantamount  to  a  textbook.  This  would  include  a  discussion  of  the 
physi.. a  involved  and  the  various  means  for  determining  the  numerical 
values  of  each  of  the  individual  hydrodynamic  coefficients  required  for 
the  equations.  In  addition,  it  is  planned  to  discuss  the  various  computer 
techniques  used  for  simulation  and  to  present  those  computer  programs 
which  nave  reached  an  advanced  state  of  development.  Such  a  treatise, 
ho  we « er,  is  considered  to  be  much  beyond  the  scope  of  what  is  required 
to  satisfy  the  immediate  needs  stated  in  the  request. 

To  accomplish  the  specific  objectives  mentioned  at  the  outset,  this 
report  presents  a  brief  history  of  the  development  of  the  equations  of 
motion  at  NSRDC;  defines  the  nomenclature,  axes  systems,  and  sign 
conventions  used  for  the  equations;  presents  the  NSRDC  Standard 
Equ  auons  of  Motion  (or  conducting  submarine  simulation  studies;  outlines 
the  methods  and  sources  for  obtaining  the  hydrodynamic  coefficients  and 
other  input  data  required  for  simulation;  and  briefly  discusses  the 
application  and  range  of  validity  of  the  Standard  Equations. 

HISTORY 

The  equations  of  motion  now  in  use  at  NSRDC  are  the  direct  result 
of  a  continuing  program  of  ln-house  research  carried  out  by  the  Stability 
and  Control  Division  of  the  Hydromechanics  Laboratory.  This  program 
was  officially  started  in  June  1956  when  it  became  evident  that  a  rapid 
technological  expansion  was  taking  place  in  the  area  of  submarine 
stability  and  control.  The  decision  was  made  at  that  time  to  develop 
general  equations  of  motion  for  six  degrees  of  freedom  which  could  bo 
used  in  conjunction  with  coefficients  obtained  from  model  experiments 
to  predict  and  study  the  behavior  of  submarines  and  associated  systems 


in  a  variety  of  maneuver*.  Concurrently,  studies  were  undertaken  to 
develop  suitable  techniques  to  determine  experimentally  all  of  the  hydro - 
dynamic  coefficients  required  for  these  equations  for  any  arbitrary  sub¬ 
merged  body -appendage  configuration.  The  Planar -Motion -Mechanism 
System,  which  was  constructed  and  placed  into  regular  service  in  19S7, 
was  one  outgrowth  of  these  studies.  This  system  is  still  the  standard 
method  used  at  NSRDC  to  determine  numerically  most  of  the  required 
coefficients.  It  was  realized  at  the  outset  that  many  years  of  research 
would  be  required  to  develop,  implement,  and  validate  a  completely 
general  set  of  equations  that  could  be  used  to  study  all  situations  involved 
in  the  motions  of  submerged  bodies.  In  view  of  the  pressing  needs,  how¬ 
ever,  it  was  decided  to  approach  the  ultimate  objective  by  progressive 
stage#. 

The  first  effort  was  directed  toward  the  development  of  equations  to 
permit  simulation  of  normal  types  of  maneuvers  involving  six  degrees  of 
freedom,  such  as  those  required  by  submarines  to  fulfill  effectively  their 
various  missions.  The  existing  linearized  equations  of  motion  were  first 
used  as  a  basis  to  study  moderate  or  small-scale  maneuver e.  These 
were  later  extended  to  include  various  nonlinear  and  coupling  term#.  By 
1960,  it  was  felt  that  the  state  of  the  art  had  progreascd  to  the  point  where 
#t  ne  standardization  would  be  desirable.  Accordingly,  the  StabJUty  and 
Control  Division  of  the  Hydromechanics  Laboratory  adopted  a  "standard" 
set  of  tix-degrce-of-f  reedom  equations  for  submarines  in  normal  types 
of  maneuvers.  The  term  standard  was  applied  in  the  sense  that  the 
equations  were  to  be  used  for  their  designated  purpose  in  all  studies 
made  by  the  Laboratory,  and  were  to  have  official  status  until  a  suf¬ 
ficient  advam .e  was  made  in  the  state  of  the  art  to  justify  a  change.  If 
such  turned  out  to  be  the  case,  it  was  intended  that  the  revisions  would 
be  made  officially  and  a  new  set  of  standard  equations  would  be  issued. 

Together  witv  the  adoption  of  standard  equations,  a  continuing  program 
v*a*  initiated  lo  correlate  full-scale  trial  measurements  with  computer 
piedictions  involving  these  equations.  On  basis  of  correlation  studies 
made  subsequently,  it  appeared  that  these  standard  equations  together 
with  coefficient#  determined  with  the  Planar -Motion -Mechanism  System 
wruld,  tor  most  part,  *,  '.eld  accurate  predictions  of  the  submarine 
trajectories  for  a  variety  o.  normal  or  definitive  maneuvers  thi.  could 
be  used  to  evaluate  the  handling  qualities  of  submarines. 

During  1962,  in  anticipation  of  a  program  directed  toward  the 
development  of  control  proceoures  and  auxiliary  device©  to  improve 
the  emergency  reco-er >  capabilities  of  submarines,  the  research  program 
was  reoriented  to  establ  '*>  prediction  techniques  fer  studying  the  behavior 
of  submarines  in  a  variety  of  possible  emergency  situations.  The  revised 
research  program  included  'undamental  theoretical  and  exper'  \ental 
studies  such  as:  expanding  the  standard  equations  of  motion  Include 
term#  (s  .ch  as  those  associated  with  "backing"  on  propellers  and  blowing 
of  bellast  tanks)  which  were  not  required  for  simulation  of  normal 
maneuvers;  extending  the  range  of  validity  of  hydrodynamic  cocfiicienta 
to  provide  accurate  representation*  in  extreme  situations  (ouch  a*  those 
involving  angles  of  attack  up  to  90  degrees);  and  est«bu#hing  correlation 
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between  predicted  trajectories  and  those  obtained  on  full-scale  submarines 
during  "emergency  recovery  trials". 

About  December  1963,  the  Bureau  of  Ships  requested  the  Laboratory 
to  consider  undertaking  a  program  of  model  tests  and  simulation  studies 
to  determine  the  ability  of  each  of  twenty-two  different  submarine  types  to 
recover  from  sternplane  jam  emergencies.  By  this  time,  considerable 
progress  pertinent  to  the  sternplane  jam  problem  had  been  made  on  the 
development  of  the  equations  of  motion.  Therefore,  the  Laboratory 
agreed  to  accept  full  responsibility  for  carrying  out  the  proposed  program, 
provided  the  simulation  studies  could  be  carried  out  with  equations  that 
constituted  the  "ctate  of  the  art"  as  of  February  1964. 

In  July  1964,  at  the  request  of  the  Bureau,  the  Laboratory  extended 
its  program  to  include  load  supportability  studies  on  the  various  submarine 
types.  This  study  was  well  within  the  scope  of  existing  equations  and 
consequently  no  further  development  was  required.  However,  in  Septem¬ 
ber  1964,  when  the  problem  of  roll  associated  with  free  rise  of  submarines 
after  blowing  was  introduced,  it  became  necessary  to  develop  new  model- 
test  techniques  and  to  conduct  further  analytical  studies  to  extend  the 
equations  to  permit  detailed  simulation  studies  of  the  problem.  Meanwhile, 
as  the  result  of  research  studies  on  the  Rotating  Arm  Facility,  further 
refinements  were  made  in  certain  coupling  terms  associated  with  "squat¬ 
ting"  behavior  while  turning,  and  these  were  incorporated  into  the  equations. 
Thus,  the  equations  contained  In  this  report  represent  the  state  of  the  art 
which  has  existed  since  about  March  1965. 


STANDARD  EQUATIONS  OF  MOTION 

The  NSRDC  Standard  Equations  of  Motion  for  Submarine  Simulation, 
briefly  designated  herein  as  the  Standard  Equations,  are  presented  In 
their  entirety  in  this  section  to  provide  a  means  of  ready  reference.  The 
Standard  Equations  are  sufficiently  comprehensive  to  define  the  trajecto¬ 
ries  and  responses,  in  six  degrees  of  freedom,  of  a  submarine  undergoing 
even  extreme  maneuvers  such  as  an  emergency  recovery  from  a  hard-dive 
sternplane  jam,  including  the  buoyant  accent  after  blowing  of  ballast. 
Ordinarily,  the  Standard  Equations  are  used  for  ahead  motion;  however, 
by  pertinent  changes  in  the  coefficients  they  can  be  used  for  separate 
studies  of  astern  motions. 

A  complete  notation  for  the  various  terms  in  the  equations  is  given  in 
the  front  of  the  report.  This  notation  is  generally  in  accordance  with  the 
Standard  Nomenclature  given  in  References  2  and  3.  The  Standard 
Equations  are  referred  to  a  right-hand  orthogonal  system  of  moving  axes, 
fixed  in  the  body,  with  its  origin  located  at  the  center  of  mass  CG  of  the 
body.  The  xx  plane  is  the  principal  plane  of  symmetry  (vertical  center - 
plane  for  submarines);  the  x  axis  is  parallel  to  the  baseline  of  the  body. 

The  positive  directions  of  the  axes  are  specified  as  follows:  x-forward, 
y-starboara,  and  r-downward.  The  remaining  sign  conventions  follow 
from  the  right -hand- screw  rule.  The  positive  directions  of  the  axes, 
angles,  linear  velocity  components,  angular  velocity  componecto,  forces, 
and  moments  are  shown  by  arrows  In  the  sketch  of  Figure  1. 
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Sketch  Showing  Positive  Directions  of  Axes,  Angles,  Velocities,  Forces,  and  Moments 
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IDENTITY  AND  SOURCE  OF  INPUT  DATA 


As  mentioned  in  the  Introduction,  it  is  beyond  the  scope  of  this  report 
to  go  into  %  detailed  treatment  on  the  use  of  computer  techniques  to  perform 
simulation  studies  involving  the  Standard  Equations.  The  purpose  of  this 
section  is  merely  to  identify  the  kinds  of  input  data  required  to  perform 
such  studies  and  to  Indicate  where  such  data  can  be  obtained.  In  addition, 
an  attempt  is  made  to  indicate  some  of  the  eub- routines  that  are  also 
needed  to  perform  various  types  of  simulation  studies.  In  particular,  those 
sub -routines  required  for  the  sternplane  jam  phase  of  the  Submarine 
Safety  Program  are  Listed.  A  so/newhat  more  detailed  account  of  the 
nature  of  the  hydrodynamic  coefficients  and  the  processes  for  obtaining 
them  is  relegated  to  the  next  section. 

The  NSRDC  Standard  Equations  of  Motion  given  in  the  preceding 
section  are  written  in  a  form  utilizing  nondimen  atonal  coefficients,  and 
are  applicable  in  a  general  sense  to  the  rigid  body  motions  of  submarines 
and  other  submerged  vehicles.  Certain  Input  data  are  required,  however, 
to  obtain  a  mathematical  model  which  can  be  used  to  realistically 
simulate  the  motions  of  a  specific  submarine  design.  A  complete  set  of 
input  data  for  this  purpose  consists  of  the  numerical  values  of  all  of  the 
hydrodynamic  coefficients,  inertial  properties,  and  pertinent  geometric 
characteristics  for  the  given  submarine  that  enter  into  the  various  terms 
of  the  equations.  Fortunately,  as  the  result  of  the  concentrated  effort  on 
the  Submarine  Safety  Program.as  well  as  contract  design  studies  made 
prior  to  construction,  complete  sets  of  inp  data  for  about  25  different 
existing  submarine  designs  are  on  hand  at  N..ADC.  The  data,  for  most 
part,  r»rry  the  classification  of  "CONFIDENT!  AD'.  However  upon  request, 
they  c  .n  be  furnished  to  outside  agencies  and  contractors  having  the 
proper  security  clearance  and  "need  to  know4’. 

Similar  data  can  be  obtained  for  new  submarine  designs  by  address¬ 
ing  a  request  to  the  Commanding  Officer  and  Director  of  the  Naval  Ship 
Research  and  Development  Center.  This  request  should  be  accompanied 
by  an  outline  of  program  objectives,  lines  plans,  and  other  pertinent  data. 
An  estimate  of  time  and  cost  will  be  provided,  after  which  funds  must  be 
deposited  at  NSRDC  before  work  can  proceed.  If  the  information  is  needed 
on  a  short  term  basis,  this  should  be  clearly  stated  by  the  requesting 
agency  so  that  sufficient  lead  time,  particularly  to  allow  for  model  con¬ 
struction  and  scheduling  of  test  facilities,  can  De  allowed.  Furthermore, 

In  contractual  dealings,  an  eariy  commitment  to  supply  the  data  should 
be  obtained  from  NSRDC  before  final  negotiations  are  made.  This  is 
especially  stressed  for  those  contracts  in  which  time  is  of  the  essence, 
or  where  the  contractor  may  understand  that  having  prior  possession  of 
the  information  is  a  condition  precedent  to  performance. 

The  primary  mathematical  model,  consisting  of  the  Standard 
Equations  and  input  data  pertinent  to  a  specific  submarine,  can  be  readily 
programmed  on  either  an  analog  or  digital  computer.  To  perform  a 
complete  motion  simulation  study,  bowsver,  the  computer  program  must 
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also  include  certain  supplementary  mathematical  model*  or  sub-routines. 
For  example,  to  conduct  »  complete  simulation  study  of  the  motions  in 
six  degree*  of  freedom  of  a  submarine  involved  in  emergency  recovery 
from  a  siernplane  jam  casualty  require*  at  least  the  sub-routine*  which 
include  the  following  representations: 

1.  Time  history  of  sternpl&ne  movement  in  degrees  including  various 
prescribed  sternplane  "jams" 

Z.  Time  history  of  sailplane  movement  in  degrees  starting  from  time 
the  order  is  given 

3.  Time  history  of  rudder  movement  in  degrees  starting  from  time 
the  order  is  given  with  provision  for  various  modes  of  motion 
such  as  swinging  the  rudder  sharply  over  to  a  hardover  position 
(either  right  or  left),  or  fishtailing 

4.  Time  history  of  change  in  propeller  rpm  from  "full  ahead"  to 
"back  emergency"  starting  from  the  time  the  order  is  given  and 
including  appropriate  time  allowances  for  communication,  human 
reaction  and  manipulation  of  controls,  and  prime -mover  response 

5.  Time  history  of  weight  of  ballast  wat^r  dlscarged  or  taken  on  due 
to  blowing  or  venting,  respectively  each  of  the  main  ballast 
tanks.  The  time  is  measured  from  the  instant  the  order  is  given 
and  includes  the  effects  of  system  lags,  as  well  as  pressure 
change*5  due  to  changes  in  air -bank  pressure  and  submarine 
depth  on  rate  of  water  discharge. 

The  first  three  of  these  sub-routines  convert  time  histories  of 
control  surface  movement  to  time  histories  of  hydrodynamic  force*  and 
moments  exerted  on  the  submarine  hull  when  they  are  applied  to  the 
terms  in  the  primary  mathematical  model  containing  8,,  6r,  and  6». 

The  time  history  of  propeller  rpm  is  converted  to  external  hydrodynamic 
forces  and  moment*  by  means  of  the  terms  containing  (rj-1)  and  u„. 

These  terms  reflect  the  incremental  changes  in  forces  and  momenta  due 
to  either  over  or  under  propulsion  (including  backing  on  the  propellers 
white  the  submarine  is  proceeding  ahead).  For  the  moderate  changes 
in  ahead  speed  involved  in  most  normal  maneuvers;,  all  of  the  (17- 1 ) 
terms  usually  can  be  neglected.  The  sub-routine  for  blowing  or  venting 
of  ballast  tanks  is  already  In  the  form  of  forces  and  moments.  It  can  be 
tied  in  directly  with  the  primary  mathematical  model  by  replacing  the 
constant  weight  terms,  such  as  VV-B,  with  the  appropriate  time  histories 
of  weight  changes. 

It  Is  apparent  that  the  primary  mathematical  model,  In  conjunction 
with  the  foregoing  types  of  sub -routines,  can  be  used  to  simulate  the 
submarine  undergoing  a  host  of  other  normal  submerged  maneuvers  as 
well.  These  maneuvers,  however,  are  essentially  of  the  open  loop  type. 
To  perform  closed  loop  studies,  such  ae  depthkeeping  or  coursekeeping 
under  the  influence  of  various  environmental  conditionj,  other  sub-routines 
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including  either  manual  operators  or  automatic  control  systems  rttust  be 
added.  A  typical  simulation  study  of  this  type  is  the  problem  of  main¬ 
taining  periscope  depth  while  operating  under  a  heavy  sea  (say  State  5). 

A  sub-routine  to  represent  the  effects  of  an  ahead,  random  sea  on  the 
submarine  has  been  developed  and  used  by  NSRDC  on  a  number  of  oc¬ 
casions  in  the  past.  This  sub-routine,  sometimes  called  the  sea-state 
analog,  is  in  the  form  of  forcing  functions  which  are  generally  applied  in 
an  analog  computer  simulation  by  means  of  a  white-noise  generator,  the 
output  of  which,  in  some  cases,  has  been  prerecorded  on  magnetic  tape. 

The  major  components  of  the  sea-state  analog  are  as  follows: 

1.  The  surface  configuration  of  the  representative  random  seaway 

2.  The  effects  of  the  seaway  on  the  depth  (pressure)  sensor  which 
are  functions  of  depth 

3.  The  oscillatory  forces  and  moments  acting  on  the  submarine  hull 
which  are  functions  of  depth 

4.  The  suction  forces  acting  on  the  submarine  which  are  functions  of 
depth 

The  Standard  Equations  have  been  programmed  at  NSRDC  both  on  the 
analog  computers  of  the  Systems  Simulation  Facility  of  the  Hydrome¬ 
chanics  Laboratory  and  the  digital  computer  (IBM  7090)  of  the  Applied 
Mathematics  Laboratory.  These  programs,  however,  have  not  yet  been 
adequately  documented  for  general  publication.  Similarly,  computer 
programs  exist  for  a  variety  of  sub-routines.  Since  these  sub-routines 
usually  must  be  tailored  to  a  particular  problem,  or  in  some  cases  to  a 
particular  computer,  no  concerted  effort  has  been  made  to  systematize 
and  document  them  for  general  publication.  Nevertheless,  the  input  data 
required  to  construct  the  sub-routines  for  existing  submarines  which  have 
been  studied  are  available  at  NSRDC  and  can  be  furnished  to  those  entitled 
to  receive  them.  Other  sources  for  such  input  data,  particularly  those 
for  sub- routines  dealing  with  hardware  items  such  as  propulsion  machinery, 
ballast  blowing  systems,  control  linkages,  and  automatic  control  systems 
are:  the  Naval  Ship  Engineering  Center  (NAVSEC),  the  submarine  con¬ 
struction  shipyards,  and  the  equipment  manufacturers. 


THE  HYDRODYNAMIC  COEFFICIENTS 


The  hydrodynamic  coefficients  constitute  the  heart  of  the  primary 
mathematical  model  used  in  simulation  studies  of  rigid-body  motions  of 
submarines.  In  fact,  the  form  selected  for  the  equations  of  motions  is 
influenced  to  a  large  extent  by  the  kiod3  of  hydrodynamic  coefficients 
employed.  With  this  qualification  in  mind,  this  section  of  the  report 
presents  further  information  as  to  the  nature  of  the  hydrodynamic  coef¬ 
ficients  used  at  NSRDC  it  well  as  some  insight  into  bow  they  are  obtained. 

It  is  hoped  that  this  will  serve  to  assure  that  the  hydrodynamic  coefficient* 
being  used  for  a  patticv'lar  study  are  compatible  with  the  Standard  Equations. 
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Unless  mentioned  otherwise,  the  hydrodynamic  coefficient*  apply  to  the 
deeply  submerged  case,  free  of  free-surface,  bottom,  and  wall  effects. 

The  hydrodynamic  forces  and  moments  which  enter  into  the  Standard 
Equations  as  coefficients  are  usually  classified  into  three  general  cate¬ 
gories:  static,  rotary,  and  acceleration.  The  static  coefficients  are  due 
to  the  components  of  linear  velocity  of  the  body  relative  to  the  fluid;  the  j 
rotary  coefficients  are  due  to  components  of  angular  velocity;  and  the  ‘ 

acceleration  coefficients  are  due  to  either  linear  or  angular  acceleration 
component*.  Within  limited  ranges,  the  coefficients  are  linear  with 
respect  to  the  appropriate  variables,  and  thus  may  be  utilized  as  static, 
rotary,  and  Acceleration  derivatives  in  linearized  equations  of  motion. 

For  six  degree*  of  freedom,  the  coefficients  in  each  of  the  three 
categories  which  appear  in  the  Standard  Equations  are  quite  numerous. 

The  problem  thus  becomes  one  of  determining  the  numerical  values  of 
the  individual  coefficients  with  sufficient  accuracy  to  support  the  objectives 
of  the  particular  simulation  study  involved.  Ideally,  it  would  be  desirable 
to  acquire  the  required  numerical  values  for  a  given  submarine  configura¬ 
tion  by  means  of  hydrodynamic  theory.  Unfortunately,  those  coefficients 
which  are  primarily  due  to  viscous  flow,  such  as  the  static  and  rotary 
coefficients,  cannot  be  obtained  reliably  using  only  existing  theory.  Theory 
has  been  used  with  reasonable  success  to  compute  acceleration  coefficient 
for  simple  forms  without  appendages,  such  as  bodies  of  revolution,  whidy 
are  amenable  to  treatment  on  basis  of  potential  flow  considerations.  How¬ 
ever,  for  actual  submarine  configurations  which  include  appendages  j«ch 
as  control  surfaces,  decks,  bridge  fairwaters,  and  propellers,  thcu&se  of 
theory  to  determine  the  values  for  coefficients  even  of  this  type  wnfi  the 
accuracy  required  for  certain  simulation  studies  becomes  somewhat 
questionable.  Accordingly,  the  present  state  of  the  art  is  to  jmy  almost 
entirely  on  experimental  means  for  determining  the  numerioftt  values  of 
the  hydrodyrv.  mic  coefficients  for  specific  design  problem^/ 

The  primary  methodology  used  at  NSRDC  to  obtain  i*fe  required  hydro¬ 
dynamic  coefficients  is  the  Flanar-Motion-Mechanisn^Systein3.  This 
systc  m  was  conceived  and  developed  in  the  StabilittOTnd  Control  Division 
of  the  Hydromechanics  Laboratory.  It  incorporaUpif  in  one  device  a  means 
for  experimentally  determining  all  of  the  types  c/nydrodynamic  coef¬ 
ficients  in  each  of  the  three  categories  that  are^requlred  in  equations  of 
motion  for  a  submerged  body  in  six  degree*  fijy  freedom.  A  complete 
description  of  the  system  including  its  basU^concepts,  principles  of 
operation,  apparatus,  instrumentation,  typical  test  results,  is  given 
in  References  5  and  6.  Planar-Motion^Sechanism  System,  Mark  I,  as 
it  is  presently  constituted,  can  be  me/ to  obtain  the  numerical  values  of 
all  of  the  coefficients  except  certain/coupling  terms  and  those  non- 
linearities  associated  with  high  v abaci  of  nondimensional  angular  velocity 
component, /such  as  r',  which  ar^asoociated  with  maneuvers  Involving 
relatively  yght  turns .  Where  each  values  are  required,  the  data  obtained 
from  Planar-Motion-MechanUiTn  tests  are  supplemented  by  the  results  of 
Routing-Arm  tests  and/or  estimates  based  on  theory.  The  Rotating  Arm 
Facility  afnd  the  instrumentation  and  technique#  used  to  conduct  submarine 


model  tests  thereon  are  descri  ^ed  in  References  7  and  6,  respectively. 
Further  detail*  on  the  experirr  ntal  and  e*timation  technique*  which 
relate  to  the  kind*  and  quality  f  coefficient*  to  be  u*ed  specifically  with 
the  Standard  Equation*  are  presented  in  the  following  paragraphs. 

The  Standard  Equation*  ar  written  in  term*  of  the  complete  sub- 
marine  configuration.  Accord  ngly,  for  quantitative  studies  pertaining 
to  specific  submarine  designs,  the  standard  practice  at  NSRDC  is  to  con¬ 
duct  all  of  the  required  tests  v.  th  models  that  are  fully  equipped  with  all 
significant  appendages  including  bridge  fairwater,  deck,  bowplanes  or 
sailplanes,  siernplanes,  rudders  and  propellers.  These  tests  effectively 
cover  a  matrix  containing  a  range  of  motion  variables  (referred  to  the  hull 
center  of  gravity)  that  generall  /  exceeds  that  which  could  be  encountered 
by  the  prototype  submarine,  ar  d  ranges  of  control  surface  angles  and  pro¬ 
peller  rpm's  which  are  at  leaf  equal  to  the  capacity  of  the  specified  sub¬ 
marine.  The  resulting  charac’eriaation  of  hydrodynamic  forces  and 
moment*  thus  embraces  interaction  effects  involved  in  the  various  modes 
of  rigid  body  motion  between:  control  surfaces  and  hull,  propeller  and 
hull,  and  propeller  and  stern  control  surfaces.  Also  included  are  down- 
wash  effects  of  forward  appendages,  such  as  bridge  fairwater,  sailplanes, 
and  deck,  acting  on  the  stern  control  surfaces. 

Since  the  foregoing  process  yields  hydrodynamic  coefficients  which 
pertain  to  the  complete  configuration,  there  is  no  need  to  include  separate 
terms  in  the  equations  in  an  attempt  to  account  for  local  angles  of  attack 
in  way  of  control  surfaces.  The  latter  procedure  has  been  suggested  by 
some  other  investigators  as  a  means  for  obtaining  a  better  representation 
of  control-surface  stall,  particularly  where  large  local  angles  of  attvck 
occur  when  the  control  surfaces  are  deflected.  On  the  other  hand,  the 
procedure  for  testing  the  complete  configuration  embraces  the  effects  on 
control  coefficients  due  to  changes  in  local  angles  of  attack  resulting 
from  the  various  modes  of  motion.  Since  true  flow  angles  of  attack  are 
inherent  in  this  process,  the  resulting  control  coefficients  should  provide 
a  more  accurate  representation.  In  fact,  if  the  two  procedures  are  com¬ 
bined,  and  numerical  values  are  assigned  to  the  various  coefficients,  there 
is  always  a  possibility  of  creating  a  redundancy  in  the  equations.  It  is 
strongly  recommended,  therefore,  if  a  complete  set  of  hydrodynamic 
coefficients,  based  primarily  on  experiment  is  furnished  by  NSRDC,  that 
the  numerical  values  of  all  coefficients  be  used  directly  without  alteration. 

In  addition  to  their  being  completely  equipped  and  self-propelled,  as 
a  matter  of  standard  practice,  the  models  used  at  NSRDC  for  tests  involv¬ 
ing  specific  submarine  designs  are  large  (usually  about  20  feet  in  length). 
Such  large  model*  used  in  conjunction  with  the  large,  rigid  towing  tank 
facilities  of  the  Hydromechanics  Laboratory  permit  the  determination 
of  hydrodynamic  coefficients  which  are  corriparatively  free  of  scale 
effects  and  other  extraneous  experimental  problems.  The  advantages  of 
large  models  and  large  facilities  are  particularly  prevalent  when  obtain¬ 
ing  data  at  high  angles  of  attack  on  the  hull  and/or  control  surfaces.  It  is 
essential,  in  both  c  as  re,  that  the  Reynoldp  number  U  t  be  high  enough  to 
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Where  fairly  complete  data  are  available  for  a  specific  submarine 
dep-.^n  as  the  result  of  an  extensive  series  of  tests  with  the  Planar -Motion- 
Mechanism  System,  the  decision  frequently  is  made  to  forego  supple¬ 
mentary  Rotating- Arm  tests  to  obtain  the  missing  nonlinearities  mentioned 
previously.  This  is  because  estimation  techniques,  based  on  empirical 
data  derived  from  prior  Rotating-Arm  tests  with  comparable  submerged 
bodies,  have  been  developed  by  NSRDC.  Indications  to  date  are  that  these 
estimation  techniques  are  sufficiently  accurate  for  the  purpose  of  making 
nearly  all  of  the  various  kinds  of  simulation  studies  required  for  most 
modern  types  of  military  submarines.  Therefore;  in  these  cases,  it  is 
felt  that  any  possible  further  refinement  would  not  justify  the  costs  of  the 

supplements ry  *.er*g  with  the  Rotating-Arm  Facility.  Also,  as  mentioned 
previously,  there  are  a  few  coupling  coefficients  that  are  not  being 
presently  obtained  by  either  of  the  two  experimental  techniques.  These 
are  estimated  on  basis  of  theory  using  the  relationships  given  in  Reference 
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avoid  the  effects  of  transitional  flow  over  the  hull  and  appendages  of  the 
model.  The  effects  of  transitional  flow  at  large  hull  angles  of  attack  are 
manifested  as  cross -flow  drag  coefficients  that  are  much  higher  than  they 
should  be  for  the  corresponding  full-seale  submarine.  The  effect  of 
transitional  flow  on  the  control  surfaces  is  primarily  manifested  as  prema¬ 
ture  stall  (breakdown  in  forces)  on  the  model  control  surfaces.  This  results 
in  control  coefficients  at  the  larger  control  deflections  which  are  lower 
than  they  should  be  for  the  corresponding  full-scale  submarine.  With  the 
large  models,  a  sufficiently  high  Reynolds  number  can  be  obtained  at  a 
moderate  speed,  say  about  6  knots  for  a  20-foot  model,  to  avoid  scale 
effects  in  most  cases.  With  the  large  facilities,  i.e.  facilities  having  cross- 
sectional  dimensions  that  are  large  compared  with  the  corresponding 
dimension  of  the  model,  these  Reynolds  numbers  can  be  obtained  without 
encountering  blockage  or  free-surface  effects  which  require  corrections 
to  the  data.  Other  advantages  in  the  use  of  large  models  include:  minimi¬ 
zation  of  tow-strut  interference  effects,6  reduction  of  propeller  scale 
effects,  ease  of  model  alignment,  accuracy  of  model  ballasting,  and  ability 
to  house  standard  instrumentation,  propulsion  motors,  and  control  actuators 
within  the  model.  All  of  these  factors  contribute  to  the  accuracy  and  repeat¬ 
ability  of  the  end  results. 

it  should  be  noted  that  the  use  of  small  models  (5  feet  long  or  less) 
does  not  in  itself  preclude  the  ability  to  obtain  hydrodynamic  force  and 
moment  data  which,  for  the  deep  submergence  case,  are  essentially 
free  of  transitional  flow  effects.  For  example,  in  water  of  the  same 
temperature,  the  same  Reynolds  number  can  be  obtained  with  a  5-foot 
model  at  a  speed  of  24  knots  as  that  mentioned  as  adequate  for  the  20- 
foo;  model  at  6  knots.  Unfortunately,  most  of  the  laboratories  which 
employ  small  models  have  relatively  small  facilities  especially  if  they 
have  the  speed  capability.  Therefore,  the  blockage  effects  constitute  the 
limiting  problem,  particularly  when  conducting  tests  at  large  hull  angles 
of  attack. 
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The  hydrodynamic  cata  required  for  simulation  studies  of  a  specific 
submarine  design  are  usually  transmitted  by  NSRDC  to  the  users  in  a 
standardized  tabular  form.  The  table  contains  a  complete  listing  of 
the  numerical  values  for  each  of  the  coefficients  shown  in  the  Standard 
Equations.  The  numerical  values  of  the  coefficients  are  derived  as 
follows: 

The  hydrodynamic  force  and  moment  data  obtained  by  the  afore¬ 
mentioned  experimental  and  analytical  techniques  are  first  reduced  to 
nondimensional  form  in  accordance  with  the  normalization  formulas 
given  in  the  Notation.  These  data  can  then  be  plotted  or  tabulated  as 
functions  of  each  of  the  appropriate  nondimensional  kinematic  variables. 

In  general,  a  least  square  fit  can  then  be  applied  to  the  data  representing 
each  of  the  functional  relationships  to  obtain  the  desired  nondimensional 
coefficients.  Certain  functional  relationships  are  known  to  be  linear  from 
considerations  of  theory,  such  as  those  associated  with  the  acceleration  or 
"added  mats"  coefficients.  In  such  cases,  only  a  straightline  fit  is  made, 
and  the  resulting  numerical  values  can  be  used  interchangeably  as  coef¬ 
ficients  in  the  Standard  Equations  or  as  stability  derivatives  in  linearized 
equations  of  motion.  For  nonlinear  functional  relationships,  least  square 
fits  to  polynomials  are  used.  The  current  practice,  in  nearly  all  cases, 
is  to  carry  such  fits  only  up  to  the  second  order  coefficients.  This  is 
done  primarily  to  facilitate  computer  representation,  particularly  where 
analog  computer  techniques  are  to  be  employed.  The  use  of  a  second 
order  representation,  however,  has  at  least  some  foundation  in  theory 
for  certain  quantities  such  as  those  arising  from  cross-flow  drag.  It 
should  be  emphasized  that  the  numerical  values  of  first  order  coefficients 
associated  with  second  order  fits  are  not  necessarily  equal  to  those 
customarily  used  for  corresponding  stability  and  control  derivatives  in 
linearized  equations  of  motion.  However,  the  standard  notation  used  to 
describe  the  various  first  order  coefficients,  including  the  acceleration 
coefficients  previously  mentioned,  is  the  same  as  that  used  for  the  cor¬ 
responding  stability  and  control  derivatives.3  To  avoid  misunderstanding 
therefore,  all  of  the  values  contained  in  the  standardized  tablet  provided 
by  NSRDC  are  to  be  taken  as  coefficients  of  the  Standard  Equations  to  be 
used  for  simulation  studies,  unless  mentioned  otherwise.  Any  table  which 
contains  numerical  values  being  provided  ior  use  in  stability  and  control 
analyses  involving  linearized  equations  oFmotion,  will  be  clearly  marked 
as  a  table  of  stability  and  control  derivatives. 

/ 

It  may  be  noted  in  the  Standard  Equations  that  some  of  the  coefficients 
involve  absolute  values  of  certain  kinematic  variables.  This  is  done  to 
assure  the  proper  signs  in  the  computer  representations  and,  in  some 
cases,  to  obtain  a  better  representation  when  a  function  happens  to  be 
asymmetrical,  such  as  the  variation  of  M*  with  w*  with  respect  to  rise 

and  dive.  / 

/ 

Although  the  standardized  tables  list  all  of  the  coefficients  shown  in  the 
Standard  Equation*,  some  of  the  coefficients  may  be  assigned  a  value  of 
zero.  This  is  because,  lor  some  submarine  types,  the  values  are  actually 
zero  or  small  enough  to  be  neglected.  Furthermore,  for  certain  submarine 
typos,  if  may  not  be  necessary  to  represent  the  coefficients  that  vary  with 
T)  when  studying  only  normal  ♦vpes  of  maneuvers. 
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In  addition  to  the  hydrodynamic  coefficients  shown  explicitly  in  the 
Standard  Equations  for  the  deep  submergence  case,  NSRDC  provides 
other  hydrodynamic  data  for  uae  in  submarine  simulation  sb-dies. 

Typical  data  in  this  respect  are  the  coefficients  or  forcing  functions 
associated  with  proximity  to:  the  free  surface  (either  in  still  water  or 
waves),  the  ocean  bottom,  or  otltcr  boundaries.  These  data  are  obtained 
by  model  experiment,  by  theory,  or  by  a  combination  of  doth. 

An  example  of  such  supplementary  data  are  the  hydrodynamic  forcing 
functions  used  to  construct  the  sea-state  analog  which  was  mentioned 
earlier  as  one  of  the  sub-routines  used  for  depthkeeping  studies.  These 
forcing  functions  are  compounded  on  the  basis  of  the  forces  and  moment# 
acting  on  th  submerged  body  moving  under  trains  of  regular  waves  of 
progressively  varied  frequencies  and  amplitudes.  When  a  body  moves  at 
constant  speed  beneath  a  regular  train  of  waves  of  a  given  period  and 
amplitude,  two  general  types  of  forces  arise:  oscillatory  forces  having  a 
frequency  equivalent  to  the  encounter  frequency,  and  a  constant  force 
called  a  suction  force  which  always  tends  to  pull  the  submarine  toward 
the  surface.  The  oscillatory  forces  and  moments  are  determined  ex¬ 
perimentally,  in  some  cases.  This  is  done  by  conducting  special  tests 
in  which  the  model  is  restrained  to  the  towing  carriage  at  each  of  several 
depths  near  the  surface  and  towed  over  a  range  of  constant  speeds  for  each 
of  si  vtral  uniform  wave  conditions.  The  resulting  forces  and  moments 
and  the  wave  configuration  are  recorded  as  time  histories.  The  phase  angle 
between  the  forces  and  moments  and  wave,  referred  to  the  hull  center  of 
gravity,  is  also  determined.  Oscillatory  forces  and  moments  which,  in 
manv  cases,  are  sufficiently  accurate  for  depthkeeping  studies  can  be 
determined  by  the  method  of  Reference  9.  The  suction  forces  are  dif¬ 
ficult  to  determine  experimentally.  Therefore,  the  usual  practice  is 
to  compute  them  by  meins  of  theory. 10 

Once  the  oscillatory  and  suction  forces  have  been  determined  for  a 
range  of  cases  involving  regular  waves,  they  are  compounded  into  spectra 
for  the  case  of  a  representative  random  sea,  say  a  State  5  sea,  by  a  technique 
developed  by  the  Hydromechanics  Laboratory.  It  should  be  mentioned  that 
although  the  suction  forces  generated  by  this  technique  are  random,  they 
always  remain  directed  toward  the  surface.  The  data  for  the  seaway  con¬ 
figuration,  the  forcing  functions,  and  response  of  the  depth  sensor  to  the 
seaway  are  presented  in  spectral  form,  including  functions  to  account 
for  attenuation  with  increase  in  depth.  In  some  cases,  the  data  have  been 
recorded  on  magnetic  tapes  for  runs  of  20-minute  duration  which  can  be 
used  directly  in  an  analog  computer  simulation. 


APPLICATION  AND  RANGE  OF  VALIDITY 


The  NSRDC  Standard  Equations  of  Motion  can  be  used  in  conjunction 
with  the  aforementioned  hydrodynamic  coefficients  and  sub-routines  to 
perform  a  wide  variety  of  submarine  simulation  studies.  Several  of  these 
applications  have  been  mentioned  earlier  in  a  somewhat  different  context. 
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The  purpose  of  this  eection  is  to  present  a  more  complete  listing  of  the 
various  applications  and  discuss  their  ranges  of  validity. 

Simulation  techniques,  such  as  those  discussed  in  this  report,  have 
progressed  to  the  point  where  they  have  become  the  primary  tools  for 
studying  rigid-body  motions  and  related  phenomena  pertaining  to  sub¬ 
marines.  Aa  mentioned  in  the  Introduction,  such  techniques  already 
have  been  applied  to  solve  a  wide  variety  of  problems  pertaining  to  the 
design  and  operation  of  submarines.  In  addition,  they  have  been  used 
effectively  with  training  simulators  to  achieve  improved  Fleet  readiness, 
with  large  savings  in  personnel-training  costs.  Undoubtedly,  there  are 
numerous  other  applications  which  will  become  apparent  as  time  goes  by. 
The  following  is  a  compilation.,  not  necessarily  all  inclusive,  of  the  various 
categories  of  simulation  studies  that  have  been  carried  out  using  the 
Standard  Equations.  Unless  mentioned  otherwise,  all  of  the  maneuvers 
involved  in  these  studies  per'iin  to  the  submerged  condition  in  ahead 
motion. 

1.  Definitive  maneuvers  {open  loop)  to  evaluate  inherent  handling 
qualities 

a.  Meanders  (vertical  plane) 

b.  Vertical  overshoots 

c.  Horizontal  steady  turns 

d.  Horizontal  overshoots 

e.  Horizontal  spirals 

f.  Acceleration  and  deceleration  in  straightline  motion 

2.  Normal  maneuvers  (operational  or  tactical) 

a.  Depthkeeping  and  coursekeeping  at  various  speeds,  in¬ 
cluding  hovering,  using  manual  or  automatic  control 

1.  Deeply  submerged  with  environmental  disturbances 
such  as  density  gradients  and  cross-currents 

2.  Near-surface  under  various  sea  states 

3.  Near-bottom  including  large-scale  bottom 
irregularities 

b.  Limit  dives  using  manual,  semi-automatic,  or  automatic 
control 

c.  Transient  horizontal  turns  using  manual,  semi-automatic, 
or  automatic  control 

d.  Spiral  descents 

e.  Mission  profiles  of  various  types  including  target 
tracking,  weapons  delivery,  and  a  variety  of  evasive 
maneuvers 

3.  Emergency  Maneuvers 

a.  Recovery  from  sternplane  jam  casualties  using  various 
combinations  of  recovery  measures 
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b.  Recovery  from  flooding  casualties  using  various  combinations 
of  recovery  measures 

c.  Buoyant  ascents  to  develop  safe  procedures  for  exercising 
emergency  ballast  blow  systems 

d.  Maneuvers  to  determine  load  supportability  as  a  function  of 
speed 

The  main  strength  of  the  simulation  techniques  described  herein 
lies  in  their  ability  to  predict  or  represent,  with  reasonable  accuracy, 
the  behavior  of  a  given  submarine  in  a  variety  of  modes  of  motion. 
Furthermore,  if  maximum  benefit  is  to  be  derived,  these  techniques 
should  perform  this  function  either  in  advance  of  construction,  or  with¬ 
out  involvement  of  the  full-acale  submarine  where  it  already  exists.  To 
arrive  at  this  state  of  the  art,  it  is  apparent  that  an  active  and  vigilant 
correlation  program  to  verify  the  accuracy  and  validity  of  predictions 
must  be  maintained. 

The  Hydromechanics  Laboratory  devotes  a  large  part  of  its  efforts 
to  the  business  of  prediction,  namely  providing  solutions  to  problems 
b^'rre  the  fact,  and  therefore,  has  recognized  the  need  for  a  stronp 
correlation  program.  Such  a  program  has  been  maintained  on  a  continu¬ 
ing  basis  by  the  Stability  and  Control  Division  over  the  past  ten  years. 

The  overall  objective  of  this  program  is  to  determine  the  accuracy  with 
which  the  stability  and  control  characteristics  of  submarines,  surface 
ships,  and  other  marine  vehicles  can  be  predicted  by  various  alternative 
model  and  analytical  techniques.  In  the  more  recent  years,  a  large  seg¬ 
ment  of  the  program  has  been  devoted  to  establishing  correlation  between 
computer  predictions  based  on  techniques  of  the  type  described  herein, 
and  measurements  taken  on  full-scale  trials  of  submarines. 

To  provide  meaningful  correlation  data  for  submarines,  the  need 
for  designing  and  carrying  out  carefully  conducted  special  types  of 
full-scale  trials  cannot  be  overemphasized,  in  addition,  to  assure  succbbs, 
these  special  trials  should  be  carried  out  as  a  complete  package  and  not 
intermingled  with  operational  maneuvers  or  other  runs  being  performed 
for  entirely  different  objectives.  It  is  beyond  the  scope  of  this  report  to 
go  into  detail  as  to  the  manner  in  which  these  trials  are  conducted.  Suf¬ 
fice  it  to  say  that  the  subject  submarine  must  be  highly  instrumented  with 
accurate  sensors  and  recording  equipment  which  may  either  supplement 
or  be  used  in  lieu  of  some  of  the  ship's  own  instrumentation.  Sufficient 
instrumentation  coverage  must  be  provided  not  only  to  measure  the  direct 
kinematic  quantities  involved,  but  also  some  of  the  more  subtle  quantities, 
which  may  affect  later  comparisons.  Furthermore,  in  conducting  such 
trials,  strict  attention  must  he  paid  to  controlling  initial  or  reference 
conditions  including  such  precautionary  measures  as  periodically  taking 
"stop  trims". 

For  a  number  of  years,  the  Stability  and  Control  Division  has  been 
conducting  full-scale  trials  primarily  to  evaluate  and  establish  numerical 
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measures  of  handling  qualities  of  submarines.  In  general,  these  trials 
consist  of  definitive  maneuvers  of  the  types  listed  previously.  The 
definitive  maneuvers  are  programmed  maneuvers  carried  out  under 
controlled  conditions  in  which  measurements  of  all  significant  motion 
variables  are  taken  with  accurate  instrumentation.  Consequently,  they 
are  ideal  for  providing  correlation  data  which  are  representative  not 
only  of  definitive  maneuvers,  but  most  normal  operational  maneuvers  as 
well.  Included  in  the  handling  quality  trials  are  some  depthkeeping  and 
coursekeeping  runs  at  periscope  depth  under  various  sea  states.  The 
data  from  these  tests  can,  at  best,  be  used  as  a  statistical  tie-in  with 
the  co.Tputer  predictions.  Since  the  foregoing  types  of  trials  are  usually 
carried  out  on  the  first  of  each  new  class  of  submarines,  correlation 
data  representative  of  normal  maneuvers  have  been  acquired  on  most  of 
the  different  types  of  submarines  that  have  been  put  into  service  from  the 
ACSS  5 69  on. 

Correlation  data  related  to  the  various  modes  of  motion  associated 
with  emcr  -  ,  cy  recovery  maneuvers  have  been  considerably  less  abundant. 

This  is  because  the  emergency  recovery  maneuvers  which  were  included 
in  the  handling-quality  trials  were  conducted  primarily  to  demonstrate  the 
submarine's  capability  of  recovering  from  sternplant  jams  of  varying 
degrees.  Consequently,  accurate  measurements  were  not  taken  of  the 
various  events  and  other  quantities  pertinent  to  the  correlation  problem. 
Accordingly,  about  years  ago,  it  became  necessary  to  design  and 
initiate  a  series  of  special  trials  specifically  directed  toward  providing 
the  desired  correlation  data.  These  trials  were  conducted  under  the 
sponsorship  of  the  Submarine  Safety  Program.  They  were  of  two  types: 
one,  to  provide  data  on  the  behavior  of  representative  submarines  in 
sternplane  jam  recovery  maneuvers;  and  the  other,  to  provide  data  on 
behavior  during  buoyant  ascents  following  blow  ing  of  the  main  ballast 
tanks.  Obviously,  in  both  types,  the  program  was  confined  to  maneuvers 
that  could  be  carried  out  without  endangering  the  submarine,  as  pre¬ 
determined  on  the  basis  of  computer  predictions.  Up  to  the  present  time, 
correlation  trials  of  both  the  sternplane  jam  and  emergency  blow  variety 
have  been  conducted  on  about  six  different  submarines,  although  not 
necessarily  the  same  ones  In  each  case. 

It  will  be  several  yearo  before  all  of  the  data  accumulated  from  the 
special  trials  of  the  various  submarine  types  can  be  thoroughly  analysed 
and  issued  in  a  formal  report  containing  a  comprehensive  treatment  of 
the  correlation  problem  for  all  of  the  modes  of  motion  involved.  In  the 
interim,  it  is  planned  to  issue  separate  reports  for  individual  submarines 
which  will  compare  measured  trajectories  with  computed  trajectories, 
based  on  the  Standard  Equations,  for  both  normal  and  emergency  maneuvers. 
These  reports  should  serve  to  provide  an  indication  of  the  current  status 
of  the  prediction  techniques.  Some  reports  of  this  type  have  been  issued 
for  the  cases  of  normal  vertical-plane  maneuvers,  r~rmal  horizontal- 
plane  maneuvers,  and  buoyant  ascents  associated  with  emergency  blows. 

In  addition,  comparisons  of  measured  and  predicted  trajectories  were 
made,  preparatory  to  conducting  simulation  studies  to  determine  safe 
operating  limits  and  emergency  recovery  capabilities,  for  about  ten  dif¬ 
ferent  submarine  types.  These  comparisons  included  representative 
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normal  maneuver*  auch  as  vertical  overshoots,  horizontal  turns,  and 
deceleration  runs  as  well  as  some  moderate  sternplane  jam  emergency 
recovery  maneuvers. 

On  the  basis  of  the  correlation  studies  and  preliminary  comparisons 
made  to  date,  it  appears  that  the  Standard  Equations  together  with  coef¬ 
ficients  of  the  type  described  herein  will,  for  most  part,  yield  accurate 
predictions  of  all  pertinent  trajectories  associated  with  a  variety  of 
normal  definitive  or  operational  maneuvers  in  submerged  ahead  motion. 

In  fact,  the  agreement  in  the  majority  of  the  cases  investigated  was  well 
within  the  ability  of  the  full-scale  submarine  to  repeat  the  same  maneuver. 
For  emergency  recovery  maneuvers,  generally  good  agreement  has  been 
found  between  measured  and  computed  trajectories  for  moderate  stern- 
plane  jam  recoveries  and  buoyant  ascents  after  bal'aat  blows.  Some  un¬ 
usual  effects  and  discrepancies  have  been  discovered  while  conducting 
emergency  recovery  maneuvers  on  a  few  of  the  submarines.  Most  of  the 
discrepancies  have  been  traceable  to  a  variety  of  factors  other  than  the 
hydrodynamic  coefficients.  However,  certain  unusual  effects  such  as 
gross  differences  in  behavior  between  the  use  of  right  and  left  rudder 
as  a  recovery  measure,  have  not  yet  been  adequately  explained.  In  sum¬ 
mary,  within  the  range  of  the  correlation  studies  made  to  date  and  in 
spite  of  the  few  discrepancies  noted,  there  appears  to  be  no  valid  reason 
to  ct.ange  the  Standard  Equations  as  they  now  exist. 

The  complete  range  of  validity  of  the  Standard  Equations  cannot  be 
established  solely  on  the  basis  of  correlation  studies  of  the  types  described 
herein.  At  best,  these  studies  are  limited  to  maneuvers  that  can  be  under¬ 
taken  by  the  submarine  with  a  reasonable  margin  of  safety.  For  example, 
ccrmal  maneuvers  for  high-speed  submarines  are  not  carried  out  to  their 
lull  potentiality.  Also,  the  emergency  recovery  maneuvers  are  generally 
carried  out  at  much  less  extreme  conditions  than  those  that  could  be  en- 
cou-iicred  in  real  casualties.  Consequently,  it  is  necessary  to  use  a 
somewhat  different  basis  to  assess  the  validity  of  the  Standard  Equations 
beyo  d  the  range  covered  by  the  correlation  trials. 

Baa’cally,  the  range  of  validity  of  the  Standard  Equations  could 
extend  to  all  modes  of  motion  encountered  by  a  submarine  during  sub¬ 
merged  fli'ht  »t  zero  or  any  ahead  speed.  This  would  embrace  various 
normal  and  emergency  situations  involving  angles  of  attack  or  drift  up  to 
90  degrees.  Tc.  achieve  this  range,  however,  would  require  a  complete 
characterisation  of  the  hydrodynamic  forces  and  moments  acting  on  a 
submarine  model  wh’ch  covers  a  broad  matrix  of  kinematic  variables  and 
combinations  thereof,  including  various  control-surface  deflections.  The 
characterizations  currently  being  made,  particularly  those  being  used  for 
the  Submarine  Safety  Program,  are  fairly  extensive  and  include  angles  of 
attack  and  drift  up  to  90  degre-s.  The  results  of  such  characterisations 
are.  for  most  part,  incorporate  in  the  hydrodynamic  coefficients  in  the 
Standard  Equations.  However,  a  simplified  representation  is  still  being 
uned  for  the  control  coefficients  associated  with  the  deflection  of  stern- 
planes,  bowplanes  or  sailplanes,  anu  rudder.  Although  simplified,  the 
representation  for  the  control  coefficient;?  appears  to  give  approximately 


24 


the  same  result*  a*  are  obtained  by  a  more  rophisticated  representation, 
even  for  fairly  extreme  maneuver*. 

The  simplified  representation  could  become  inadequate  only  if  the 
effective  angle  of  attach  at  a  control  surface  become#  large,  causing  a 
condition  of  stall  in  which  the  lift  no  longer  increases  with  the  deflection 
angle.  In  normal  submerged  maneuvers,  including  steep  dives  and  tight 
turns,  the  effective  angles  of  attack  at  the  control  surfaces  of  a  submarine 
tend  to  be  substantially  smaller  than  the  nominal  angles,  and  consequently 
stall  will  not  be  experienced.  There  are,  however,  some  unusual  con¬ 
ditions,  such  as  may  occur  in  certain  flooding  casualties,  which  could 
produce  effective  anglea  at  the  diving  planes  large  enough  to  cause  stall. 

If  the  condition  Is  transient,  and  the  large  angle  oi  attack  is  only  of  short 
duration,  there  may  be  little  effect  on  the  trajectories  by  using  the  more 
sophisticated  representation.  However,  in  studies  of  lead  supportability, 
large  angles  of  attack  can  be  obtained  on  the  sternplanes,  particularly 
where  the  load  is  concentrated  at  the  stern.  Since  thee 3  studies  involve 
a  steady-state  maneuver,  the  use  of  the  simplified  representation  for 
sternplanes  could  be  erroneous  for  the  type  of  loading  mentioned.  Rather 
than  encumber  the  Standard  Equation  to  treat  these  few  special  cases,  a 
separate  sub-routine  to  represent  diving  plane  stall  has  been  developed 
and  is  available  at  NSRDC.  The  sub-routine  represents  the  leveling  off 
of  lift,  which  is  an  approximate  characteristic  of  control-surface  stall, 
but  allows  the  drag  to  increase  with  deflection  angle. 

In  addition  to  the  unusual,  and  as  yet  unexplained,  effects  which  have 
been  observed  on  full-scale  Submarines,  there  are  several  other  factors 
which  may  be  pertinent  to  the  range  of  validity  of  the  Standard  Equations. 
Included  among  these  are  the  effects  of  vortex  shedding  on  rolling 
behavior,  and  effects  due  to  the  submarine  moving  close  to  the  free 
surface  while  carrying  large  pitch,  roll,  or  yaw  angles  of  orientation. 

A  research  program  to  investigate  these  factors  and  which  could  possibly 
lead  to  future  refinements  and  extensions  of  the  Standard  Equations,  is 
being  maintained  at  NSRDC. 


CONCLUDING  REMARKS 

The  Standard  Equations  presented  in  this  report  have  been  used 
successfully  by  NSRDC  for  the  past  Z$  years  to  conduct  simulation 
studies  of  specific  submarines  engaged  in  submerged  maneuvers  ranging 
from  normal  maneuvers  to  extreme  maneuvers  such  as  those  encountered 
in  emergency  recovery  from  sternplane  jam  and  flooding  casualties. 
During  this  period,  several  hundred  of  such  studies  involving  about  25 
different  submarine  designs  have  been  carried  out  and  the  results  of  many 
of  these  studies  are  contained  in  the  classified  literature.  A  companion 
program  of  correlation  studies  is  also  being  carried  out  to  determine  the 
extent  to  which  computer  predictions,  involving  the  use  of  the  -Standard 
Equations,  agree  with  measured  trajectories  obtained  from  full-scale 
trials  of  submarines  undergoing  various  type*  of  maneuvers.  Those 
comparisons  made  to  date  have  generally  shown  good  agreement  between 


predicted  and  measured  trajectories  of  representative  maneuvers  of 
both  the  normal  and  emergency  types.  It  can  be  expected,  therefore, 
that  the  NSRDC  Equations  used  in  conjunction  with  the  specified  hydro¬ 
dynamic  coefficients  and  tub -routines  will  provide  simulations  of 
trofectoriea  that  are  well  within  the  accuracy  required  for  the  purposes 
of  those  kinds  of'studiee  conducted  up  to  the  present  time. 

Although  the  NSRDC  Standard  Equations  have  reached  a  fairly 
advanced  stage  of  development,  they  still  must  be  considered  only  as 
repi  e tentative  of  the  state  of  the  art,  and  subject  to  change  if  so 
dictated  by  future  research.  An  active  research  program  is  being 
maintained  at  NSRDC  to  examine  some  of  the  remaining  problem  areas 
and  to  continue  the  development,  if  required.  However,  in  interest  of 
consistency,  a  policy  has  been  established  not  to  change  the  Standard 
Equations  unless  it  can  be  demonstrated  that  significant  improvements 
in  predictions  can  be  obtained.  Such  a  change  muat  be  officially  made 
by  NSRIlC,  and  all  parties  concerned  will  be  properly  notified. 

In  view  of  the  foregoing,  it  is  recommended  that  the  NSRDC  Standard 
Equations  and  associated  methodology  be  adopted  aa  a  U.S.  Navy-wide 
standard  for  simulation  of  specific  submarine  deaigna,  and  that  they  be 
given  official  status  in  dealing  with  contractors  and  other  outside 
activities. 
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